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We demonstrate the approach towards a Gibbs-like equilibrium state, with a common temper- 
ature and a chemical potential, of two finite metallic grains, prepared with a different number of 
noninteracting electrons, connected by a weak link that is susceptible to incoherent and inelastic 
processes. By developing an analytic method and by using an exact numerical approach, the quan- 
tum time evolution of the electrons in the metallic grains is followed. In the absence of decoherring 
and inelastic effects, equilibration is never reached. Introducing dephasing effects on the link only, 
using a dephasing probe, the two quantum systems equilibrate, but do not evolve towards a Gibbs- 
like state. In contrast, by mimicking inelastic interactions with a voltage probe, the metal pieces 
evolve towards a common Gibbs-like equilibrium state, with the probe. 

PACS numbers: 05.30.-d, 03.65.Aa, 03.65.Yz,72.10.-d 



Introduction. How do quantum systems equilibrate 
from a certain non-equilibrium initial condition, e.g., a 
quench [l|? With cold atoms in optical lattices offering 
a clean experimental setup , renewed attention in this 
fundamental topic has recently sparked. One could ad- 
dress this question with (at least) three distinct setups 
in mind: (i) Consider an isolated quantum system, and 
study its evolution towards equilibrium, for example, in 
the time averaged sense 0, H| . (ii) Attach a subsystem 
with a few degrees freedom to a thermal reservoir, and 
monitor the system equilibration, e.g., in the sense of the 
small trace distance or (iii) put in contact two iden- 

tical finite quantum systems and watch for the process 
of mutual equilibration Q. Despite intense efforts, gen- 
eral results are still missing P . Recent studies also argue 
about the precise definition of quantum integrability, and 
its implication on quantum thermalization [3, [13] • 

In this work, we consider a combined setup and demon- 
strate the process of mutual equilibration of two finite 
metallic quantum systems, connected by a weak link, in 
the presence of either elastic-decohcrring processes or in- 
clastic effects, through the interaction of the link elec- 
trons with additional degrees of freedom. When only 
decoherence effects are allowed, the system approaches 
a non-canonical equilibrium state. In contrast, when in- 
clastic processes are included, the two parts relax towards 
a common Gibbs-like state. The origin of inelastic scat- 
tering processes are many-body interactions in the sys- 
tem, e.g., coupling electrons to phonons. Since an explicit 
and exact inclusion of such effects is challenging pdTjl4| , 
it was suggested [H, [l6| to phenomenologically intro- 
duce elastic and inelastic scattering processes by using 
dephasing or voltage probes, respectively. These electron 
reservoirs are prepared such that there is no net (energy 
resolved or total) electron flow from the system towards 
these probes. 

As a particular realization, we consider the non- 



interacting Anderson model with a single electronic level 
(dot) coupled to two metallic grains (henceforth referred 
to as reservoirs) [l7| . Each reservoir is initially prepared 
in a distinct Gibbs-like grand canonical state, at a dif- 
ferent chemical potential. We follow the time evolution 
of the reservoirs' electrons, once put in contact through 
the dot, itself susceptible to decoherring and/or inelastic 
processes. For a schematic representation, see Fig[TJ We 
refer to the metal grains, including (each) N ~ 500 elec- 
tronic states and n ~ 200 electrons as "reservoirs", to 
indicate that they have a dense-enough density of states, 
such that their effect on the impurity (dot) can be ab- 
sorbed into a positive real self-energy function, allowing 
for a quantum Langevin equation (QLE) description fl8| . 
While we may consider large reservoirs, the number of 
electrons in the metal grains-dot unit is fixed. However, 
energy is conserved only in the elastic scattering scenario. 
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FIG. 1: Two metallic grains (reservoirs) separately prepared 
in a grand canonical - diagonal state; the initial population 
is plotted at the boundaries. At to the reservoirs are put 
into contact through a single electronic state, susceptible to 
decoherring and relaxation effects. This is accounted for by 
the G reservoir, serving as a dephasing or a voltage probe. 



Model. We consider two electron reservoirs v = L,R, 
with identical density of states and a sharp cutoff at ±D, 
which do not directly couple, only through their weak hy- 
bridization with a single-state quantum dot. The Hamil- 
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tonian takes the form 

H =H l +H r + H w + V l + Vr, (1) 

where Hl,r,w represents the Hamiltonian for left reser- 
voir, right reservoir and dot, respectively. The term V v 
denotes the coupling of the dot to the v reservoir, 

H L = y^Qqc;, H R =) j e r 4cy, H w = e d c\c d 

l r 

Vl = ^ VlC^Cl + h.C. Vr = V r c d C r + h -C. (2) 

I r 

Here, ct (ck) are fermionic creation (annihilation) oper- 
ators of the left reservoir, I 6 L, right reservoir, r £ R 
or the dot (d). We assume that vi and w r are real num- 
bers and that the Hamiltonians of the metal grains are 
diagonal in momentum basis. A factorized initial state is 
assumed, with an empty dot and the reservoirs prepared 
in a diagonal (grand canonical) state at a chemical po- 
tential p v and inverse temperature j3 = 1/T, satisfying 
the distribution f v (e) = [e /3(e ~ Al " ) + l]" 1 , f i L = -p, R . 

At t t ne two reservoirs are put into contact through 
the dot, and their dynamics is followed using either an 
exact quantum evolution scheme, or a QLE approach (de- 
tails below). The observed population dynamics is de- 
picted in Fig. EJa). When the dot level is placed within 
the bias window, a resonance pattern shows around e d 
with a peak developing in the accepting R reservoir and 
a corresponding dip showing at the L side. The dynam- 
ics shown in Fig. [Ha) is fully coherent. Results pre- 
sented are before the recurrence time r° ec ~ 2ir/AE; 
AE = 2D/N is the mean spacing between energy levels 
Q. At this time, a complete depletion of certain levels 
occurs, and the dynamics is reversed. 

We now wish to allow for elastic dephasing effects or 
inelastic interactions on the dot only. We mimic such 
effects with "probes" , by augmenting the Hamiltonian 
([2|) by an additional (finite size) electron reservoir G, 

H = H + H G +V G . (3) 

Here H G = J2 g e g c l c g and V G = J2 g Ugcjd+h.c. Inelas- 
tic effects are introduced using a voltage probe [HI, , 
demanding that the net current from the dot to the G 
unit vanishes, i G = 0. This condition sets an effective 
chemical potential for this reservoir. Alternatively, elas- 
tic dephasing effects can be introduced using a dephasing 
probe, requiring that ic(e) = 0, i.e., the charge current 
at a given energy should vanish. Results are displayed 
in Fig. HJb) and (c): The system approaches equilibrium 
(dephasing probe), and even thermal equilibrium (volt- 
age probe). These results are discussed in details below 
Eq. (fTU|) . We now explain how we time-evolved the sys- 
tem's density matrix under Hq or H. 

QLE Method. We aim at calculating the time evolu- 
tion of all two-body operators in the system. We begin 
with the trivial part, the impurity (dot). Since it is cou- 
pled to many degrees of freedom, its dynamics can be 



described using a quantum Langevin equation (l8l - [2l| . 
We review the steps involved, so as to highlight the un- 
derling approximations, setting the limit for the method 
applicability. We outline the derivation in the absence 
of the G probe; We generalize it later to include such a 
device. In the Hcisenberg representation the operators 
satisfy the following equations of motion (EOM), 

Cd = -itdCd - l v l c l - i VrC r 
l r 

ci = -ieici - iviCd, c r = -ie r c r - iv r c d - (4) 
Using a formal integration, the I operators follow 

c,(t) = e-W-^aW-im I dre-^-^Cd^dT. 

J to 

(5) 

Similar relations hold for the R operators. We plug these 
expressions into the dot EOM and get the exact result 

c d = -ie d c d - f dT^h~ ltl(t ~ T) Cd{T) 

J to ; 

- / drJ2 «, 2 e- ie " (t - T) c d (r) - ^ L (i) - i V R (t). (6) 

Jto ,, 

Here, if{t) = J2i vie~ iei{t ~ to) ci{t ) } and similarly ?y H , 
represent "noise" . Wc now assume that the second and 
third terms reduce, each, into a dissipation term, further 
inducing an energy shift of the dot energy, absorbed into 
the definition of td- This is justified here as the metal 
grains play the role of charge and energy baths with re- 
spect to the dot. Under the markovian approximation, 
we reach the (time local) QLE 

Cd = -ie d c d - ii] L (t) -in R (t) -T(e d )c d , (7) 

withr(e) = £„ =iifl r„(e) and e.g., T L (e) =nJ2 l v?6(e- 
e;). Equation Q thus relies on two assumptions: (i) 
a positive real self-energy function r„(e) can be writ- 
ten [l8|, and (ii) the dot dynamics is slow relative to 
the reservoirs' evolution. We now use the exact equa- 
tion (|5|) and the reduced result ([7j), and derive ana- 
lytic expressions for the expectation values (c k (t)cj(t)) = 

Tr[p(t )cl(t)cj(t)}; k,j = l,r,d. Here p(t ) = Pdd>PL®PR 
is the factorized time-zero density matrix of the system. 
The trace is performed over all degrees of freedom. In 
particular, the following initial condition is assumed 

(4(*o)ca(*o)> = 0, (4(toMt )) =f L (e t ) =fi, 
(4(t )c r (to)> = f R (e r ) = fr- (8) 

with f u (e) = [e^( e_AI ") + The resolved occupation 

of the, e.g., L reservoir, is given by three contributions, 

p(c,) = (4(t)c 4 (t)) = {cl(t ) Cl (t )) 
+iv l e- i ^ t - t ^ f e i£ '( t - T )(4(r) Q (to))dr + C . C . 

Jto 

+vf f f dnd^icliny^e^-^e-^-^. 

JtQ JtQ 
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The first term accommodates the initial condition. The 
second element (denoted by F 2 ) represents first order 
reservoir-dot coupling processes. The last term (Fa) 



corroborates higher order effects, including population 
transfer from the r reservoir, 



J 



F 2 = -vffi x(t- t ) 



2T 



Fa 



r 2 



Vl 



— — <^ 2 [e 2 , - T 2 ] cos [e dl (t - to)] + 4e d Tsin [e dl (t - to)] 



vffie 



£ 2 
t dl 



y k , f 4 sin 2 [Y (*- h)] , 1 



E dk' 



r 2 



£ ■ 

k'eL.R 

1 e -(t-io)(r+ie dI ) + e -(t-t )(T+ie dk ,) _ e -i(t-t )e lk , 

{e d i - iT) ei k > 



e -2T(t-t ) + 1 _ e (t-to)(ie d i-T) _ e -(t-t )(ie dl +r) 



+ c.c 



(9) 



Here, Cjk = £j — £fe. For simplicity, the energy depen- 
dence of the dissipation terms is ignored. We can immedi- 
ately confirm that i u (t) = J2keu( c k(t) c k(t)) produces 
the standard expression for the charge current [22[. We 
have similarly derived closed expressions for all density 
matrix elements, including off-diagonal elements, e.g., 
(cj(t)ci'(t)), (I ^ I'). These lengthy expressions are not 
presented here. Since the reservoirs include many states, 
after a short time rt, Trt > 2, the dot dynamics should 
remain fixed at a quasi steady-state value, up to the re- 
currence time T° ec . The interval r f < t < r® ec is identified 
as the quasi steady-state (QSS) region 

We now incorporate the G reservoir, the probe, into 
this QLE description: The dot dynamics follows Eq. 
([7]) with an additional noise term 77°, and we re-define 
the total hybridization, T = T L + T R + T G , T G (e) = 
7r^ g v 2 (5(e — € g ). As a result, all expectation values fol- 
low a form technically identical to the T G = limit. For 
example, the population (c\(t)a(t)) obeys Eq. aug- 
mented by k' G G terms in Fa. To include inelastic scat- 
tering effects of electrons on the dot, we implement a volt- 
age probe assuming f G (e) = [e /3 ( e_A ' G - ) + The chem- 
ical potential \i G is set such that i G = ^ ^2 g (c\c g ) = 
is satisfied at all simulation time. With the motivation to 
explore situations beyond the linear response regime [2(| , 
we retrieve y, G numerically, by employing the Newton- 

(fe+l) _ (fe) . / (fch/-/ /„(fcK 



Raphson method [23 1 



/Iq is the initial guess, i' G is the first derivative with re- 
spect to \xq. In principle, one should adjust \i G through- 
out the simulation, to eliminate population leakage from 
the L-dot-i? system. However, we have practically found 
that within the allowed simulation time (details below) 
we could safely assume a QSS limit. The energy resolved 
charge current into G can then be written as 



2r G T R [f G (e) - f R (e)]+T L [f G (e) - f L (e)} 
(e-e d ) 2 +r 2 



with the total current i G — J i G (e)de. Quasi-elastic 
scattering effects arc implemented within a dephasing 
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FIG. 2: (a) Coherent population dynamics of the reservoirs' 
electrons. Plotted are the L (three top lines) and R (three 
bottom lines) occupations as a function of electron energy, for 
r G = at t = (o) t = 750 (full) and t = 1500 (dashed). 

(b) Approaching a non-canonical equilibrium state with a de- 
phasing probe, L-bath (five top lines) and i?-bath (five bottom 
lines) occupation functions, (c) Approaching thermal equilib- 
rium with a voltage probe, L-bath (five right-most lines) and 
i?-bath (five left-most lines) occupation functions. In (b) and 

(c) T G = 0.4 and t = (o), t = 750 (full), t = 1500 (dashed), 
t = 7500 (dashed-dotted) and t = 15x 10 3 (dotted). The 
last two lines were generated by restarting the QLE dynam- 
ics when approaching t — m,Tr ec , m is an integer, /ig has 
been further updated at each restart point to eliminate leak- 
age. In all panels /3 = 200, T L = F R = 0.025. e d = 0, 
fiL = — /j.r — 0.2, D = 1, N = 500 electronic states at each 
reservoir, L, R, G. 



probe, by demanding that i G (e) — 0, to yield f G (e) = 

r R f R (e)+r L f L (e) 

Exact method. Even with the inclusion of the probe, 
the QLE description is still technically limited by the re- 
currence time t < T® ec oc Nl.r {Nl,r is the number of 
electronic states in a single metallic grain), though no 
actual population recurrence does show in the dynam- 
ics. This is because the validity of Eq. ([7]) is limited by 
this time, beyond which inter-reservoir recurrences, that 
may not show in the overall behavior, take place. Us- 
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ing an exact, expensive, brute force calculation, we can 
numerically simulate (A = 



CjCfcj 
iHt 



(A(t))=Tr B [p(t y m A 
d 



-iHti 



limx^o-^Tr B [pLPRPGPde lHt e XA e~ 



iHU 



(10) 



Pv = 

Such 



using the fermionic trace formula [24|. Here, 
g-PiHu-fiuNv) J z u is the partition function, 
a calculation perfectly agrees with QLE results, fur- 
ther confirming that beyond r° ec , where QLE description 
breaks down, the dynamics proceed towards equilibrium, 
before r rec oc J2i=L RG^i, see Fig. [3] We now show an 
approach towards equilibrium within t < T rec using the 
exact method, and t < r® ec using the QLE technique. 
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FIG. 3: Equilibration with a dephasing probe, using an exact- 
unitary time evolution scheme, Eq. (|10p . Nl = Nr = 100, 
No = 1500, other parameters are the same as in Fig. [2jb). 
The L (full) and R (dashed) populations are shown at t = : 
St : 9St, St = 500. The dotted lines mark the values reached 
before recurrence if No = 200. The inset demonstrates a 
slow-down in dynamics in approaching the equilibrium state. 



— 0.5 

Q. 









(a) 




i ' , 






- \ i 






; 1 ' 




Cl 




1.4 -0.2 0.2 

e 



-0.4 -0.2 0.2 
e 



FIG. 4: (a) Occupation of L bath at time t = (o) and 
t = 1500 for T G = and e d = (dotted), Fq = and e d = 0.1 
(dashed-dotted), T G = 0.4 and e d = (full), T G = 0.4 and 
ed = 0.1 (dashed). The latter two lines assume the voltage 
probe condition, (b) Same for the R side occupations. Other 
parameters are the same as in Fig. [2Jc). 

Results. We identify thermal equilibration in our peer- 
quantum system setup, by adjusting the conditions of 
Refs. 0, [lj|, demanding that: (i) The system should 
equilibrate, i.e., evolve towards some particular state, 
and stay close to it for almost all time. Furthermore, the 
equilibrium state should be (ii) independent of the dot 
properties-energetics and initial state, (hi) insensitive to 
the precise initial state of each reservoir, (iv) close to di- 
agonal in the energy basis of its eigen-Hamiltonian, and 
(v) a Gibbs state. 



In Fig. [2lb)-(c) we use the QLE method and fol- 
low the reservoirs' mutual equilibration process, using 
either a dephasing probe (b) or a voltage probe (c). 
We present the reservoirs' occupation at times t/75 = 
0, 10, 20, 100, 200. Beyond T r ° ec - 1500 the QLE descrip- 
tion breaks down; The data at later times has been gener- 
ated by restarting the simulation immediately before r° ec , 
using the L and R final diagonal-distribution as an ini- 
tial condition for a new run, with the G bath re-adjusted 
to respect the probe condition. This process can be re- 
peated, to reach a complete equilibration. It can be re- 
alized experimentally by fine-tweaking the voltage probe 
(introducing a dissipation mechanism into the dynamics) 
|25j . When only dephasing effects are allowed, the sys- 
tem is approaching a non-thermal state. Ultimately, the 
population of the two reservoirs should reach a two-step 
function with p(pn < e < pl) ~ 0.5, since excess elec- 
trons at the L side loose their phase memory, therefore, 
on average, half of them populate the R side in the long 
time limit. This equilibrium state is sensitive to the pre- 
cise details of the initial electron distribution, as energy 
redistribution is not allowed. In contrast, when inelas- 
tic effects are taking place, the system does approach a 
Gibbs-likc thermal state, a step function at zero temper- 
ature. 

Simulations with a dephasing probe using the exact- 
unitary method, Eq. (|10|) . up to r rec oc J2%=l hg^'i arc 
shown in Fig. [3J demonstrating a clear evolution towards 
an equilibrium state. We build a large G so as to delay 
recurrence. Results at earlier times do not depend on the 
size of G, reinforcing the observation that G acts as an 
agent in driving the L-R mutual equilibration process. 
QLE data with G bath tweaking nicely agrees with these 
results. In order to show the analogous behavior with a 
voltage probe, a dissipative mechanism should be intro- 
duced into the G bath, e.g., by building a hierarchy of its 
interactions with the L-R system. 

Fig. [4] proves that while under coherent evolution the 
resonance peak emerges around the energy ed, in the 
presence of a voltage probe with (large enough) Tq, the 
buildup of the equilibrium state systematically occurs 
around the equilibrium Fermi energy, independently of 
the link energetics. This holds even when the dot is 
placed outside the bias window (not shown). Using a 
smaller value for Tq, temporal features show, washed out 
gradually in approaching the equilibrium state. Analo- 
gous trends also take place when allowing for dephasing 
only. 

The thermal state should be diagonal in the energy 
eigenbasis of its Hamiltonian Q. In Fig. \5\ we display 
the density matrix (DM) pij> = (ejee), excluding diago- 
nals, with and without a voltage probe, using the QLE 
technique. This quantity is expected to oscillate in the 
long time limit since the Hamiltonian is not diagonal in 
the (local) I basis. We still show the results in this ba- 
sis, so as to manifest local j/-bath properties. A subtle 
source of complication is the fact that piy decays (be- 
fore T° ec ), even without a probe, due to the finite-bias 
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assumed as an initial condition. One should therefore dif- 
ferentiate between bias-induced and probe-induced deco- 
herence processes. There are three significant differences 
in the behavior of off-diagonal elements, with and with- 
out the probe: (i) The absolute value of the coherences, 
at a given time, is smaller when Tq ^ 0. (ii) The DM ap- 
proaches a diagonal form (strict diagonal values are not 
shown), (iii) When = 0, oscillations occur around e^. 
With the probe, contributions appear mainly around the 
equilibrium Fermi energy. 
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FIG. 5: Absolute values of the density matrix elements piy 
at t — 1500 (a) Tq = 0, (b) using a voltage probe, Fq = 0.4. 
Other parameters are the same as in Fig. [2{c), td = 0.1. 

Summary. We provided evidence, using analytical and 
numerical tools, showing that two finite quantum sys- 
tems, coupled through a weak link where electrons can 
loose phase memory or exchange energy, can equilibrate 



and even thcrmalizc. Dcphasing and voltage probes are 
cunning techniques, allowing to mimic memory loss and 
energy redistribution. We could follow the time evolution 
of the peer metals numerically-exactly or using an ana- 
lytic QLE scheme, where tweaking the probes with a new 
initial state drives the peer-metal system closer and closer 
to the equilibrium state. Our results are significant for 
several reasons: (i) We show that inelastic or dcphasing 
effects on a very small- yet essential- subset of the total 
system can drive the system towards a global equilib- 
rium state, where thermal equilibration is approached in 
the former case, (ii) The colloquial, restrictive, "nonde- 
generate energy gap" condition Q is not assumed in 
our study, (iii) Standard QLE treatments [l8| overlooked 
reservoirs' dynamics all-together. Here, we frame a new 
tool for studying the dynamics of a system composed 
of many degrees of freedom, by identifying a subsystem 
and resolving its dynamics, then using this information 
backward, to re-explore the evolution of all degrees of 
freedom. Future directions include the study of electron- 
electron interaction effects (26[, and considering a quan- 
tum dot chain between the two metal grains. Here, the 



charge current coherent-diffusive crossover [271 ] may re- 
flect itself in the energy reorganization of the reservoirs. 

This work has been supported by NSERC. M.K. 
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